Pomeron dominance in deeply virtual Compton scattering 
and the femto holographic image of the proton 
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The dominance of the soft pomeron in soft high energy scattering and the evolution to the deeply 
virtual regime, predicted by perturbation theory, allow us to reveal generalized parton distributions 
from HI and ZEUS measurements of deeply virtual Compton scattering. These distributions encode 
a holographic image of the proton, which will be presented. 
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We study deeply virtual Compton scattering (DVCS) 
off the proton at high energies and aim to reveal gener- 
alized parton distributions (GPDs) Q from HERA mea- 
surements • The GPDs then yield a 3D image of quark 
and gluon distributions at small longitudinal momentum 
fraction x. The analysis substantially relies on the soft 
pomeron dominance in the non-perturbative sector. Var- 
ious points of view on the pomeron are recalled before the 
DVCS predictions are derived in a simple, however, ap- 
proximated analytic form and the analysis is performed. 

In soft high energy scattering experiments, i.e., no fur- 
ther large scale is present, the energy dependence of scat- 
tering amplitudes is described within Regge phenomenol- 
ogy as an exchange of particles in the t (momentum trans- 
fer squared) channel, see Fig. ^a). Theoretical consid- 
erations, based on partial wave decomposition, analytic- 
ity, and crossing, predict that scattering amplitudes be- 
have at large center-of-mass energy squared s as s a ^ l \ 
where a(t) = a(0) + a't is the leading Regge trajectory in 
terms of the intercept a(0) and slope oi . Such a trajec- 
tory relates the particle spin J = a{m 2 ) with its mass. 
To explain experimental measurements, the trajectory 
a P (t) = 1.08 + 0.25i has been introduced 0. It is as- 
sociated with a hypothetical spin-one quasi particle, the 
so-called soft pomeron that carries the vacuum quantum 
numbers. No particles belonging to this pomeron trajec- 
tory could be established and so it remains exceptional. 

The notion pomeron is also applied in the phenomenol- 
ogy of (semi) hard reactions. The HERA measurements, 
where a 27 GeV electron/positron beam collides with a 
820 GeV proton one, lead to the discovery that the energy 
dependence of cross sections get steeper with growing 
photon virtuality Q 2 = — q\. Let us consider the proton 
structure function i^^Bj, Q 2 ), measured in deep inelas- 
tic scattering (DIS), at small XBy The reciprocal Bjorken 
scaling variable 1/xbj = 2 Pi -qi/Q 2 , where Pi is the pro- 
ton momentum, is nearly proportional to the photon- 
proton center-of-mass energy squared W 2 = (qi + Pi) 2 . 
i<2 possesses an almost flat l/^Bj dependence at low 
Q 2 < Qo ~ lGeV 2 and gets steeper with increasing Q 2 : 



F 2 (x Bi ,Q 2 ) ~ (1/arnj) 



HQ 2 ) 



(1) 



where X(Q 2 ) > rises with growing Q 2 . 

In Regge theory A = ap(0) — 1 is considered as Q 2 inde- 
pendent, which is incongruous with experiment. To cure 



this issue, a double pomeron exchange within a soft and 
hard pomeron, e.g., ap ald (0) ~ 1.4, was suggested 0. 
The (energy ordered) resummation of the gluon ladder 
in the BFKL approach results the perturbative pomeron, 
which indeed possesses such a large intercept How- 
ever, the theoretical understanding remains poor. 

Fortunately, perturbative QCD predicts the variation 
of observables with respect to a hard scale. F 2 is the ab- 
sorptive part of the forward Compton amplitude. Analo- 
gous to Fig. Olb), it factorizes in terms of quark q(x, Q 2 ) 
and gluon G(x,Q 2 ) densities. Their scale changes are 
governed by evolution equations, written for the Mcllin 
moments, e.g., G n (Q 2 ) — f Q dxx n ~ 1 G(x,Q 2 ), as 



Q 



dQ 2 \G 




(2) 



Here E ra = ^2 q n is the flavor singlet quark density and 
7„ are the anomalous dimensions, depending on spin 
n. The gluon-gluon and gluon-quark anomalous dimen- 
sions possess a 'pomeron' pole at spin n = 1. Solving the 
evolution equation yield an exponentiation of this pole. 
In x-space this converts to the so-called double log ap- 
proximation, see below. This resummation of large con- 
tributions might be effectively parameterized as in ([T]). 

The DIS measurements are well described by perturba- 
tive QCD, where the input parton densities are extracted 
from global fits. It has been suggested to generate dy- 
namically sea quarks and gluon densities by evolution 
0. Starting from a very low input scale, fine tuning is 
required to fit experimental data. In the small x region 
the evolution is mainly driven by the 'pomeron' pole of 
the anomalous dimensions and data are indeed consistent 
with the double log approximation 0. Both approaches 
were combined within a soft pomeron like input 8] , where 
only the initial scale, two normalization factors, 



S n (Qo) 



V s 



a r (Q) ' 



G„(Qo) 



N G 



a P (0) ' 



(3) 



and perhaps the running coupling have to be adjusted. 

The dynamical generation of parton densities for small 
XBj is perhaps the most predictive approach to DIS. Al- 
though the 'pomeron' pole in the evolution operator is es- 
sential, the resulting parton densities depend on the non- 
perturbative input, too. Since real Compton scattering 
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FIG. 1: The Compton scattering amplitude (a) at high energy and small t is dominated by the pomeron exchange in the 
t-channel. The lowest order perturbative contribution to DVCS at the scale Q 2 = —q\ and at a very low input scale is shown 
in (b) and (c), respectively. The integration path used in Eq. (1J and its deformation in the complex n-plane are shown in (d). 



can be described at high energy within soft pomeron ex- 
change Q , it is apparent that the non-perturbative input 
is dominated at a low scale by the soft pomeron, too. It 
arises in the expectation value of operators, sandwiched 
between proton states, as pole contribution ©. 

Let us apply this conjecture to DVCS. At leading or- 
der (LO) the factorization of the amplitude is given by 
the hand-bag diagram, depicted in Fig.^b), in terms of 
Compton form factors (CFFs) 0. At small £ = 2Bj/2 
only the flavor singlet part of the helicity conserved and 

elevant 11]. To 
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charge even CFF H(^,t, Q 2 
diagonalize the evolution operator, the complex confor- 
mal partial wave decomposition will be employed |l2l Il3| 



W = — 

1% 



c+ioo 



dn 



T(n + 3/2) 
r(3/2)r(2 + r 



7rn\ 



(4) 



|) (Q|,0)E„(Q 2 ,Q2 ) 



Here the singularities in the complex n-plane are on 
the l.h.s. of the integration path, cf. Fig. dd). Q 2 S = 
12qQq/ n f i s t ne averaged charge square fraction for n/ 
active quarks. E ra and G n are conformal GPD moments, 
which reduce in the forward case to the parton density 
ones. The evolution operator E„, a 2 x 2 matrix, is the 
same as in DIS, see, e.g., Ref. 0. By means of projec- 
tors it might be decomposed in '+' and '— ' eigenmodes: 



e«(Q 2 ,Qo) 



a s {Q) 
a s {Q ) 



- ( a s (Q) 



a s (Qo) 



,(5) 



where \^ are eigenvalues of the anomalous dimension 
matrix, appearing in Eq. @. The running coupling 



a s (Q) 



4tt 



(ll-2n / /3)In(e 2 /A|,cD) 



(0) 



is parameterized by Aqcd- 

To illuminate the physics at small XBj, let us derive 
the double log approximation for Tt s ©. The dominant 
contribution comes from the most right-lying singularity, 
see Fig.^d). If ap(t) < 1, this is an essential singularity 
appearing in the '+'-mode of the evolution operator 



A,, — 



A 



+ \++0{n-l), A 



3G 



33-2 
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For < ap(t) the next-leading singularity is expected to 
be the soft pomeron pole. To simplify the analysis, we set 
ap(0) = 1. Hence, for t = this singularity joins the lead- 
ing one and both of them should be considered in Eq. 10}. 
They can be encircled by a deformation of the integra- 
tion contour, as shown in Fig.^d). The remaining back- 
ground integral (dashed integration path c') can be safely 
neglected. Next the variable transformation n = 1 + pj 
with p = [Aln(a s (Q )/a s (Q))/ln(2/O] 1/2 and Q > Q Q 
is performed, where p is considered as small. Taylor ex- 
pansion at p = yield integrals that reduce for t = to 
the modified Bessel function I a (c) with a = 0, 1, • ■ ■ : 



J a (cr, at) 
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(8) 



where a 2 = AX\n[a s (Q )/a s (Q)} ln[2/£]. The effective 
slope a — a' F /p depends on £ and Q. Note that the 
asymptotic behavior of the integral (jSJ is a independent: 



J a (cr, at) 



1 — at 



for a 
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Now the GPD moments have to be specified. In the 
presence of a hard scale the partonic content of the 
pomeron can be resolved 01- ^ i s expected that the 
pomeron is mainly made of glue, see e.g., Ref. [l^. For 
simplicity the quark distribution, including the valence 
quarks, will be neglected in the small x region. Finally, 
the CFF reads in the double log approximation as 



H 



DL 



iirnfQ 2 s 



ph(cr, at) 



a s (Q) 
a s (Q ) 



Gifot, Ql) 



pole term 



+ o( P 2 ). 



(9) 



where A + = 1 + 20n//9(33 — 2n/). Here the pole term 
comes from the soft pomeron in the '— ' mode. It is an 
unsubstantial contribution that vanishes for t = 0. The 
real part of TC is induced by subleading terms in p and 
will not be considered here. We remark that the optical 
theorem F-i — ^smH / (2ir)\ and the Callan-Gross 




FIG. 2: In the left panel the t-dependence of the differential DVCS cross section 112H is fitted to HI data for W — 71 GeV, Q 2 = 
4GeV 2 (circles, dashed) and W = 82 GeV, Q 2 = 8 GeV 2 (squares, solid). The integrated cross section with t cut = — lGeV 2 
versus Q 2 for W = 82 [89] GeV (middle) and versus W for Q 2 — 8[9.6]GeV 2 (right) are compared with HI (circles, dashed) 
[ZEUS (squares, solid)] data 0. The statistical and systematic errors are added in quadrature. 



relation Fi = 2xBj-Pi, valid to LO accuracy, leads to the 
known double log approximation in DIS: 



F 2 (x Bh Q 2 



^i^mH DL (2xBj,i-0,2Q 2 ) (10) 
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(the factor two in scale dependence is explained by jllj)- 
The interpretation of Eq. is depicted in Fig. [He), 
where the evolution operator, i.e., the (k± ordered) re- 
summation of the gluonic ladder, is inserted in the hand- 
bag diagram. Only the moment with n = 1 is needed as 
non-perturbative input. It might be interpreted as a non- 
local gluon operator with dimension three and spin one 
that is sandwiched between the proton states. Since in 
the light-cone gauge this operator formally reduces to a 
local operator, we might argue that its matrix element is 
£ independent. In accordance with quark counting rules, 
predicting the power fall at large \t\, a simple multipol 
ansatz with unknown slope Bq is used: 



Nq 



(1 - B G t/if 



(11) 



The DVCS process is accessible via the photon lepto- 
production and interferes with the Bcthe-Heitler Brcms- 
strahlungs one. For small XBj the DVCS cross section 
dominates, the interference is negligible (after integra- 
tion over the azimuthal angle), and the Bethe-Hcitler 
cross section, known in terms of the elastic proton form 
factors, can be subtracted Q. Theoretically, the DVCS 
cross section for small x^ is safely approximated by 01 



AlTOT 



(12) 



A fit of the cross section, using Eqs. to HI 

and ZEUS data provide the parameter set 

AT G = 1.97, Bq = 3.68 GeV~ 2 , Q 2 , = 0.51 GeV 2 ,(13) 

where rif — 3 and Aqcd = 150 MeV were fixed. Note 
that the i-slope is rather robust under a change of Aqcd , 
however, it is sensitive to the form factor shape. For an 
exponential ansatz it decreases to B B q V = 2.58 GeV -2 . 



Let us confront these findings with DIS data in the 
kinematic region that is considered here. Within the pa- 
rameters (ll.il) the double log approximation l|10l) reason- 
ably describes the shape of data. However, the normal- 
ization can fail of up to 30%. Radiative QCD corrections 
are important [l(i ] and it is was stated that the normaliza- 
tion discrepancy in DVCS is resolved at next-to-leading 
order [l?|. A closer look to this problem, beyond next- 
to-leading order, will be presented somewhere else |l8|. 

The dominance of the pomeron allows us now to re- 
veal the singlet quark and gluonic GPDs, which are rep- 
resented by a Mellin-Barnes integral [12. 11^ . In the 
double log approximation they read, e.g., for vanishing 
skewness, as 



xH g (x,t] = 0) ~ I (a,at) 



a s (Q) 



xff E (x,77 = 0)~ -£ph(a,at) 
y 



a s (Qo 

a s (Q) 



Gi(t)+---, (14) 

1 A+ 



a s (Q \ 



Gi(t) 



where cr(^ = 2x) and a(£ = 2x) is set, see below Eq. |JS). 
They reduce for t — to the common parton densities. 

For 77 = the Fourier transform of H with respect to 
the transverse momentum transfer Aj_ = (-P2 — P i)±, 



rf 2 Aj 



-Ai,Q 2 ), (15) 



has a probabilistic interpretation jl9| . In the infinite mo- 
mentum frame the proton is viewed as a disc with a mean 
square charge radius r 2 = 4^-i 7 'i(i)| t _ ~ (0.6 fm) 2 . The 
density of partons is given by p(x,b, Q 2 ), depending on 
momentum fraction x, transverse distance b from the 
proton center, and resolution 1/Q. The mean square dis- 
tance for a parton species is proportional to the i-slope 



(b z ) = AB . 



B = ^- t lnH(x, V = 0,t,Q 2 ] 



t=o 



(16) 



In Fig. |3{a) the averaged distance (6 2 ) 1 / 2 is displayed 
for both gluons and sea quarks. Compared to the charge 



4 



. 95 




10"° 2 • W~ b 




0.5 



1.5 



2 

\b\ [fm] 



\b\ [fm] 



FIG. 3: The root of the mean square distance (llrJII for gluons (thick) and singlet quarks (thin) and the corresponding parton 
densities (1151 versus \b\ are shown in the left and middle panel, respectively, for x = 10~ 3 . Q 2 is set to 2 GeV 2 (dashed), 
10 GeV 2 (dash-dotted), and 100 GeV 2 (solid). The dotted lines in the middle panel display the result for an exponential 



^dependence. In the right panel the transverse distribution of giuons is visualized, where x = 10 



■ 3 and Q 2 



10 GeV 



radius r ~ 0.6 fm, this distance is of about 40-50% larger. 
Within an exponential form factor ansatz this inflation 
would shrink to 20-30%. The weakly grow with decreas- 
ing x and fall with increasing Q arise from the intertwin- 
ing of soft-pomeron input and evolution. The double log 
asymptotics, i.e., (JSJ, 114(1. and (|16|) . shed light on this: 



B(W, Q 2 ) ~B G + 



a' T In* (4W 2 /Q 2 ) 
VAln^ (a s (Q )/a s {Q)) 



(17) 



Such a behavior was also experimentally observed in hard 
vector meson electroproduction, dominated by the gluon 
exchange. Moreover, the measured i-slope of the cross 
section is consistent with our findings for DVCS |2jj. In 
Fig. I^b) it can be realized that the sea quark density is 
a bit broader and as known from forward parton densi- 
ties much smaller than the gluon one. The gluon density 
within the exponential form factor ansatz is displayed by 
dotted lines. For \b\ < 1.5 fm the density only slightly 
differs from that of the multipol ansatz. A further charac- 
teristic is that the densities (start to) vanish for \b\ > 1.5 
fm, see panel (c). However, the multipol ansatz induces 
a long tail that results into an enlargement, compared to 
the exponential one, of the averaged distance. 

In this letter we argued that the soft pomeron trajec- 
tory essentially determines the gluonic GPD at a small 
resolution scale and that perturbative evolution can be 
used to arrive at higher scales, where factorization holds 
true. This interplay of non-perturbative and perturbative 
QCD leads to a realistic behavior of the DVCS amplitude. 
To illuminate this, the double log approximation has been 
derived and used to LO. It was checked that a numeri- 
cal treatment, with changed normalization and scale pa- 
rameter, essentially leads to the same results. The fit 
to experimental data provides the non-perturbative pa- 
rameter and reveals so the 3D gluon and sea quark dis- 
tributions inside the proton. The theoretical framework 
can be refined and confronting it with precise DVCS data 
would strongly constrain the non-perturbative ansatz. 
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